X-rays, 7-rays and neutrinos from collisoinless shocks in supernova wind breakouts 



Boaz Katz 1 *, Nir Sapir 2 , and Eli Waxman 2 
^Institute for Advanced Study, Princeton, NJ 08540, USA and 
2 Dept. of Particle Phys. & Astrophys., Weizmann Institute of Science, Rehovot 76100, Israel 

We show that a collisionless shock necessarily forms during the shock breakout of a supernova 
(SN) surrounded by an optically thick wind. An intense non-thermal flash of < 1 MeV gamma rays, 
hard X-rays and multi-TeV neutrinos is produced simultaneously with and following the soft X-ray 
breakout emission, carrying similar or larger energy than the soft emission. The non-thermal flash 
is detectable by current X-ray telescopes and may be detectable out to 10's of Mpc by km-scale 
neutrino telescopes. 
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It has long been suggested that an intense burst of 
X-ray radiation is expected to be emitted at the ini- 
tial phases of a SN explosion, once the radiation me- 
diated blast wave reaches the edge of the star pH^. If 
the star is surrounded by a sufficiently optically thick 
shell of circum-stellar matter (CSM), e.g. a high density 
wind, the breakout occurs within the shell. Several ob- 
served 7-ray/X-ray flashes associated with SNe @,@| have 
been suggested [6t-TTll| to be such wind breakouts of fast, 
v sh ^ O.lc, shocks, in which departure from equilibrium 
may imply very high electron temperatures reaching tens 
or hundreds of keV [1,112] ■ In fact, all low luminosity 7- 
ray bursts associated with SNe may be produced by such 
fast, v s h O.lc, breakouts (with or without the presence 
of an optically thick CSM) HUGS- 

Breakout outbursts 
of slower shocks, u s h ~ 0.03c, have been sug gested to 
account for strong optical/UV transients |13| and very 
luminous SNe [e.g. Il4l - tla ]. In order to explain the high 
energy (reaching 10 51 erg) emitted in these SNe, CSM 
parameters were suggested such that the diffusion time 
scale through the CSM is comparable to the dynamical 
time scale, R/v s h [e.g. Ibl. Il7l - tl9| . If true, the observed 
emission is, by construction, the breakout outburst from 
the CSM. 

Following breakout, the radiation mediated shock is 
expected to become a collisionless shock, leading to the 
emission of gamma-rays and neutrinos [20| • In the ab- 
sence of a (significant) wind, the small mass of the shell 
shocked by the collisionless shock implies that only a 
small fraction, < 10 -2 , of the breakout energy is con- 
verted to such high energy radiation. Moreover, the for- 
mation of a collisionless (or collisional) shock is contro- 
versial [e.g. [3M5I [2lT - [23| . since the light shell may be accel- 
erated to sufficiently high velocity by the escaping radia- 
tion. In this letter we show that if the progenitor is sur- 
rounded by an optically thick CSM, e.g. a dense wind, a 
collsionless shock is necessarily created during the break- 
out, and that an energy comparable to or greater than 
the breakout energy is emitted by quasi-thermal particles 
in high energy (> 50 keV) photons, and by accelerated 



protons in high energy (> 1 TeV) neutrinos. The latter 
is an extension of the study of high energy emission from 
the interaction of the ejecta with a dense optically thin 
CSM 13. 

Formation of a collisionless shock. Consider first for 
simplicity a piston moving with a constant velocity v = 
10 9 cm sec -1 vg through an optically thick fully ionized 
hydrogen wind with a density profile 



p(r) = - 
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where i?b r = 10 Ru cm is a normalization parameter 
with dimensions of length. A shock propagates ahead of 
the piston with velocity i> s h ~ v. As long as the optical 
depth across the shock transition region, A r s h ~ c/v, 
is much smaller than the optical depth of the system, 
A T = (c/v)Rh x /T, the post-shock radiation is confined. 
Once the shock reaches r ~ i?br, the width of the shock 
becomes comparable to the size of the system and a sig- 
nificant fraction of the post-shock energy can be emitted 
during one dynamical time scale, Rbr/v. This emission 
is the breakout outburst discussed above. 

The material lying ahead of the piston must be accel- 
erated to velocities approacing v by some process. At 
large optical depth, where the radiation mediated shock 
is sustained, the radiation accelerates the material by 
Compton scattering off the electrons. The maximal ve- 
locity to which a fluid shell can be accelerated by this 
process is given by 



E 7 /c ctt 
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where E 7 — J Ljdt is the radiation energy emitted 
through the fluid shell and r is its initial position. This 
maximum velocity is achieved if all of the flowing photons 
move radially and the shell does not expand considerably 
during the passage of the radiation. In this case, a frac- 
tion CTr/(47rr 2 ) of the momentum E^/c carried by the 
radiation is transferred on average to each proton. 

Ery is limited by the thermal energy accumulated in the 
post shock region which, in turn, is limited by 0.5M(r)v 2 , 
where M(r) — At:(c/ Kv)Rbr r is the wind mass inward of 
r. Using equation © an upper limit for u max , w max < 



2 



0.5 (Rbr/r)v, is obtained. This implies that beyond a ra- 
dius of 0.5i?br the shock can no longer be mediated by 
radiation and must be transformed into a collisional or a 
collisionless shock. Since the ion plasma frequency, u p — 

sec -1 , is many orders 



1/2^-1/2 



(4 7 rpe 2 /m2) 1 /2 _ 10 9 w - 

of magnitudes larger than the ion Coulomb collision rate 
per particle, vc — pacv/m p ~ 2 x 10 _2 i?j" 4 1 v^" 4 sec -1 
[e.g. |2Q|, the shock will be collisionless, i.e. mediated by 
collective plasma instabilities. 

Let us comment on our simplifying assumptions of con- 
stant piston velocity and wind CSM profile. In reality, 
the shock velocity slowly changes with time as the faster- 
parts of the SN ejecta are slowed down by the wind. The 
arguments above hold if the constant v is replaced by 
its value at the vicinity of i?b r , defined as the radius at 
which t ~ c/v s h(r). Next, note that the formation of 
the collisionless shock is not restricted to the assumed 
density profile, p oc r~ 2 . Consider the more general case 
of an ejecta with mass M e j and velocity v e j propagating 
into a dense CSM with (accumulated) mass profile M{r). 



The shock propagates with velocity v s \y 



as long as 



M e j > M(r). Once the shock reaches a point where the 
optical depth is less than c/t> j breakout occurs. The con- 
siderations above imply that CSM at radii r ^> i?br for 
which M(r) -c M c j will have to be accelerated to v ~ v C j 
by a colliionless shock. 

Emission from thermal electrons. The collisionless 
shock heats the protons on a time scale of cj^ 1 to a tem- 
perature roughly given by 



Tn 
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0.2t> 2 MeV. 



(3) 



The electron temperature depends on the unknown 
amount of collisionless heating. A lower limit for the 
electron temperature can be obtained by assuming that 
there is no collisionless heating. 

The collisional heating rate of the electrons due to 
Coulomb collisions with the protons is given by 

dT p dT e 2 m e T e , 3 / 2 



The time it takes the protons to lose a significant frac- 
tion of their energy (which is of the order of the total 
available energy) is 
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The proton cooling time is thus much shorter than the 
dynamical time Rhr/v = (c/v) 2 /(naxc), where n = p/m p 
is the proton number density in the pre-shocked region 
and is smaller than rid, by the compression factor. This 
is not surprising. While the shock is radiation mediated, 
radiation energy equal to the mechanical energy is gen- 
erated on each shock crossing time scale. At breakout, 
the shock crossing time scale equals the dynamical scale 
and radiation with energy density comparable to the to- 
tal energy density must be generated during the dynam- 
ical time scale. In fact, since the electron temperature is 
higher than that expected in a corresponding radiation 
mediated shock, the emission efficiency is even higher. 

The shock is strongly radiative and the energy is effi- 
ciently converted to radiation. The typical photon ener- 
gies are expected to be of the same order of magnitude as 
the electron energies, i.e. > 60 keV. The calculation of 
the emitted spectrum is beyond the scope of this paper. 
We note that since the initial photon energies are much 
lower (~ 1 eV assuming equilibrium) we expect that the 
spectrum hardens continuously with time and that on the 
breakout time scale, significant emission is likely emitted 
at all intermediate energies. 

We conclude that gamma-rays/hard X-rays will be 
emitted with total energy comparable to that of the 
breakout energy 
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on a time scale similar to the breakout time t ~ R/v 
l-Ri4W^" 1 d with typical luminosity 



10 45 i? - a 



i4W 9 erg sec 



(9) 



where \ ep — 30A ep 1.5 is the Coulomb logarithm and 
it was assumed that T e /(m e c 2 ) ^> T p /(m p c 2 ). The 
fastest possible cooling source for thermal electrons is 
Inverse Compton (IC) scattering of the local radiation 
field, which carries a significant fraction e 7 < 1 of the 
post shock energy and is given by 



dTe 
dt 
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where C/ 7 = e-yHdTd is the photon energy density and nd 
is the shocked material proton density. Assuming U 1 
ndTp (equivalently, e 7 < 1) we find 
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When the shock expands, it will remain radiative be- 
yond i?br and the total emitted energy, integrated over 
longer time scales, may be significantly larger than that 
of the breakout energy. 

Accelerated protons: Non-thermal emission energy. 
Relativistic particles (CRs) accelerated in the collision- 
less shock that forms due to the collision of the SN ejecta 
with dense interstellar material may emit high energy 
gamma rays and neutrinos due to the interaction with 
the dense material (24|. The collisionless shock that was 
shown above to be produced during breakout from a 
dense wind is a constrained example of such interaction 
and may be a source of detectable high energy neutrinos 
and gamma rays. Here we focus on the emission from ac- 
celerated protons and their products. In what follows it 
is assumed that the accelerated protons carry a fraction 
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e CR — O.XecR,—i of the post shock energy and have a flat 
power law energy distribution, e 2 dn/d e oc e°. 

The cooling time of a relativistic accelerated proton 
due to inelastic pp collisions is roughly given by 



(Tpp \ C J V 
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(10) 



Hence, for slow enough shock velocities, v/c < 0.1, pro- 
tons accelerated at breakout efficiently convert their en- 
ergy to neutrinos, gamma-rays and pairs by pion pro- 
duction and decay (and muon decay). In this section we 
restrict the discussion to v/c < 0.1. For such shock veloc- 
ities, the amount of energy emitted by relativistic protons 
during breakout is expected to be roughly a fraction ccr 
of the energy emitted by the thermal particles. Using 
Eq. © we have ^Non-thermal ~ 10 48 e C i? ~iR 2 4 v 9 erg. 

At later stages, t pp v/r grows linearly with r, and as 
long as it is smaller than unity, the energy converted into 
pions increases linearly with the accumulated mass. The 
radius r pp at which the proton energy loss time is equal 
to the dynamical time, t pp v/r pp = 1, is 



10Wg- 2 i?br 



(11) 



Beyond this radius, the fraction of energy converted to 
pions drops like 1/r (t pp oc p -1 oc r 2 while the available 
energy increases linearly with r) implying a logarithmic 
increase in the total emitted energy. Given that in reality, 
v(r) is slowly declining, the total contribution to the non 
thermal fluence from r > r pp is of order unity compared 
to fluence produce up to this radius. The total emitted 
energy is therefore given by 
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erg. 
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Accelerated protons: Maximal proton energy. The 
maximal proton energy is limited by the time available for 
acceleration which is the shorter of the dynamical time 
and energy loss time. The acceleration time depends on 
the unknown magnetic field value and the loss time de- 
pends on the unknown target photon energy distribution. 
Nevertheless, we next demonstrate that protons are very 
likely to be accelerated to at least multi-TeV energies. 

Assuming Bohm diffusion, the acceleration time to en- 
ergy e is given by 



ta.ee — 
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where B is the post shock magnetic field and es = 
B 2 /(8irpv 2 ) is roughly the fraction of postshock en- 
ergy carried by it. For TeV CRs, the acceleration 
time is thus much shorter than the dynamical time 
and the pp energy loss time. For protons in the range 
10 — 1000 TeV the strongest possible cooling mechanism 
is photo-production of pions, with cooling time 



t pl = (0.2n 7 (jp 7 c) 1 > 5 



cr p7 m p c 2 v 
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where n 7 (/ii/ 7 ) is the target photon number density 
(typical energy) and we conservatively assumed that 
n 7 = pv 2 /(hv 1 ). Photo-production of pions occurs 
if the proton energy is higher than the threshold, ~ 
m^m p c 4 /(h^) ~ 0.13(/u/ 7 / MeV)- 1 TeV. The possible 
presence of many ~ 1 MeV photons implies that photo- 
production may be important for 1 TeV protons. Photo 
production is not important if the target photons have 
~ leV energies, as assumed in (24|. Given the constraint 
?i 7 < pv 2 1 (hi/j), the strongest losses for protons of energy 
e = £tcV TeV occurs for target photons having a typical 
energy of hv^ ~ m^m^/e ~ 0.13 MeVe^y. Using this 
in Eq. (fl4|) we obtain 
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Comparing Eq. (fl5|) to Eq. (fl3|) we conclude that ac- 
celeration to multi TeV energies is possible for e B > 
£ TcV> implying that reaching energies well 
above 1 TeV is very likely. 

We verified that proton CRs do not suffer significant 
losses due to Inverse Compton and Synchrotron emission 
during the acceleration time, and that the resulting pi- 
ons and muons do not suffer significant energy losses due 
to these processes before decaying. Finally, note that the 
maximal proton energy is increasing with radius since the 
ratio of proton acceleration time to dynamical time is in- 



dependent of radius (i acc oc B 1 



oc p 



-1/2 



oc r) while the 



ratio of all the loss times to the dynamical time decreases 
with radius. 

Accelerated protons: Multi TeV neutrinos. Roughly 
a third of the non thermal energy Eq. (fl2|) will be emit- 
ted in muon neutrinos (and anti-neutrinos) and a sig- 
nificant fraction of this energy may be emitted beyond 
TeV energies. In the neutrino energy range of one to 
hundred TeV, the effective area for muon neutrinos of 
a Cherenkov neutrino detector like IceCube is increas- 
ing linearly with energy, approximately as 10 _6 £tcV^4 5 
where A = 10 10 Aio cm 2 is the geometrical cross-section 
of the detector. The number of muons induced by one 
to hundred TeV neutrinos is therefore independent of the 
neutrino spectrum in this range and is given by 



. -g^.i-iooTcv/lO 5 erg 
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(16) 

where d is the distance to the SN and where we opti- 
mistically assumed that 1/3 of the non thermal emission, 
Eq. (fT2|) . is in multi TeV neutrinos. 

Accelerated protons: Gamma rays. High energy 
gamma-rays and pairs with energies reaching multi TeV 
energy will be generated with a comparable rate to that 
of the neutrinos. The pairs will emit further high energy 
gamma- rays by Inverse Compton interactions with the 
radiation field. Emission below ~ 1 MeV will be mixed 
with the emission from the thermal electrons. Emission 
at a photon energy hv > MeV may be suppressed by 
the large optical depth for pair creation, which depends 
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on the density of photons with energies above the pair 



production threshold hv? 



> 
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An upper limit to the optical depth for pair creation 
at a given photon energy, hv, can be obtained by us- 
ing the fact that the total energy density of photons of 
any frequency is smaller than pv 2 . Assuming that the 
energy of photons per unit logarithmic frequency does 
not exceed eo.i x 10% of pv 2 , and focusing on the radius 
10i?br at which the protons are still efficiently cooled (for 
v ~ 10 9 cmsec -1 , see Eq. pT|) ), we find 
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The emitted spectrum is suppressed by at most ~ r . 
In this 'worst case scenario', such bursts will be too faint 
to be observable by high energy (hv ~ 1 GeV) gamma- 
ray detectors such as Fermi. 

Increased optical depth due to electron-positron pairs. 
A significant fraction of the energy carried by the acceler- 
ated protons may be converted to electron positron pairs. 
The number of pairs per proton is limited by the avail- 
able energy to n e ±/n < ecR,m p v 2 / [m e c 2 ) ~ 0.2e c _R,-iWg- 
The presence of pairs cannot change the optical depth 
significantly for v/c < 0.1. For very fast shocks, v ~ c, 
the presence of pairs can potentially increase the optical 
depth considerably and affect the emitted radiation. 

Discussion. We have shown that shock breakouts in 
optically thick winds will necessarily be accompanied by 
high energy radiation from a collisionless shock that in- 
evitably forms on the time scale of the breakout outburst. 

Low luminosity GRBs associated with SNe have been 



suggested to be the outbursts associated with fast shocks 
y > 0.1c breaking out of dense optically thick winds 0- 
UM. As we have shown here, a significant fraction of the 
observed radiation, or even most of it, may be gener- 
ated by the collisionless shock that will form during the 
breakout. 

If the slow, v/c ~ 0.03, breakout interpretation of 
events such as PTF09u [13| is correct, a significant 
amount of energy, E ~ 10 51 erg, is expected to be emit- 
ted in hard X/7-rays reaching energies h v > 50 keV, 
Eq. ©, and multi-TeV neutrinos [see also[2J|. X-rays 
at lower energies are likely to be emitted with similar ef- 
ficiency and would be easily detected by instruments like 
the X-ray telescope (XRT) on board Swift or the Chan- 
dra X-ray observatory. TeV neutrinos may be detectable 
by experiments like IceCube, see Eq. (fTrJ)) , provided such 
events are sufficiently common and a similar event occurs 
at a distance d < 100 Mpc (compared to ~ 300 Mpc for 
PTF09uj). 

Finally, we note that if the CSM breakout explanation 
of very luminous SNe [e.g. [T3| is correct, our analysis im- 
plies that these events should be accompanied by strong 
high energy X-ray emission. Lack off [16], or very weak 



1 51 ] , X-ray emission from some of these events challenges 



this interpretation. 
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